Zigzag equilibrium structure in monatomic wires 
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We have applied first-principles density-functional calculations to the study of the energetics, and 
the elastic and electronic properties of monatomic wires of Au, Cu, K, and Ca in linear and a planar- 
zigzag geometries. For Cu and Au wires, the zigzag distortion is favorable even when the linear wire 
is stretched, but this is not observed for K and Ca wires. In all the cases, the equilibrium structure 
is an equilateral zigzag (bond angle of 60°). Only in the case of Au, the zigzag geometry can also 
be stabilized for an intermediate bond angle of 131°. The relationship between the bond and wire 
lengths is qualitatively different for the metallic (Au, Cu and, K) and semiconducting (Ca) wires. 



Very recently, using a combination of transmission 
electron microscopy (TEM) and scanning tunneling mi- 
croscopy, Ohnishi et al. Q were able to visualize the for- 
mation of monatomic gold wires. In particular, one of 
their images shows a wire of four gold atoms forming a 
bridge between two gold tips. The atoms in the wire were 
spaced at intervals of 3.5-4.0 A, a distance much larger 
than that of nearest-neighbors in bulk gold (2.88 A). In 
spite of these large interatomic distances, the wire was 
observed for more than two minutes, showing clearly its 
stability. Observation of gold monatomic chains with a 
length of four or more atoms was independently claimed 
by other experimental group |^] . In this work, Yanson et 
al. do not provide direct structural information, but the 
very long and flexible contacts formed in the last plateau 
of the conductance curves (with a conductance very close 
to one quantum, 2e 2 /h) were associated to the presence 
of monatomic chains. These new truly one dimensional 
structures open a great number of opportunities to test 
our understanding of the interactions in solids. However, 
the very large interatomic distances found in the experi- 
ment represent by themselves a puzzle. In ref. || we al- 
ready addressed this issue with the aid of first-principles 
density- functional calculations. Our estimate of the max- 
imum interatomic distance that the gold wires can hold 
before breaking was ~2.9 A, in clear disagreement with 
the commented experimental observations. Several re- 
cently published references also report similar conclu- 
sions |4| [?]]. The cause for such a discrepancy between 
theory and experiment was not clear. However, in ref. || 
we proposed a simple explanation for the puzzle posed by 
the TEM experiments. First, the gold wires were found 
to exhibit an equilibrium zigzag structure. Second, the 
calculated rotational energy barrier was low enough for 
the small zigzag wires, supported between gold tips, to 
be rotating at room temperature. In these conditions, 
only half of the atoms would be clearly visualized, the 
other half would produce a fuzzy image, and the appar- 
ent interatomic distance observed would correspond to 
that between second-neighbors. 



Zigzag structures for monatomic wires have already 
been observed and proposed by several authors. Whit- 
man et al. Q observed the formation of long, isolated 
zigzag chains when Cs was deposited over GaAs(llO) and 
InSb(llO) surfaces. March and Rubio Q also proposed 
zigzag monatomic chains as the building blocks for highly 
expanded liquid Rb and Cs. A zigzag distortion was also 
predicted by Chaplik |h]] for a chain of electrons held 
on the surface of liquid helium by the field of a charged 
filament. 

In this paper we use first-principles density-functional 
calculations to study the stability and electronic prop- 
erties of linear and zigzag wires, for Cu, K, and Ca, 
comparing the results with those of Au Our stan- 
dard density functional theory (DFT) flllf calculations 
have been performed with the SIESTA program |]l2| , |l3| . 
Core electrons are replaced by scalar-relativistic norm- 
conserving pseudopotentials [Q, plus a pseudocore to 
account for non-linear core corrections in the exchange- 
correlation energy jl5| . We used a doble-£ polarized basis 
set of numerical atomic orbitals, obtained from the so- 
lution of the atomic pseudopotentials with an excitation 
energy of 50 meV p6|. Real and reciprocal-space integra- 
tions were performed with a 200 Ry-cutoff grid and with 
20 inequivalent k-points, which guarantee a convergence 
better than 1 meV/atom. We present results obtained 
with the local density approximation (LDA) J!?]] for Au 
and Cu, and with the generalized gradient approximation 
(GGA) (T^| for Ca and K. However, we stress that all the 
results are very similar using LDA and GGA. We have 
performed test calculations for the bulk phases of the 
studied materials, obtaining results in good agreement 
with experiment and with previous DFT calculations. 



* The calculated equilibrium lattice parameter, bulk mod- 
ulus and, cohesive energy are: 4.11 A, 194 GPa and, 
4.40 eV/atom for Au, 3.56 A, 210 GPa and, 4.88 eV/atom 
for Cu, 5.51 A, 20.3 GPa and, 2.20 eV/atom for Ca, and 
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FIG. 1. Energy of monatomic wires of Au (a), Cu (b), Ca (c), K (d) with a planar zigzag (symbols) and a linear geometry 
(dashed lines) , as a function of the length per atom. The length x is expressed in units of the equilibrium bond length of the 
linear wires, and the energies are referred to the minimum of the linear wire. The arrows indicate points where a conduction 
band crosses the Fermi level. The triangle shows the length at which the band gap of the Ca wire moves away from the V 
point. The insests in (a) and (b) show a zoom of the energies near the equilibrium of the linear wire. 



In order to minimize the interactions between the wires 
and their periodic images, we have allowed a distance of 
20 A between wires in different cells. 

Fig [l] shows the energies of monatomic wires with a pla- 
nar zigzag and linear structure, as a function of the wire 
length, for different elements. To facilitate the compar- 
ison, lengths are given in units of the linear-chain equi- 
librium length, and energies are relative to its minimum. 
These are: 2.56 A and 2.09 eV/atom for Au, 2.25 A and 
2.16 eV/atom for Cu, 4.10 A and 0.42 eV/atom for K, 
and 4.00 A and 0.98 eV/atom for Ca. The evolution of 
the bond length in the zigzag structure can be seen in 
Fig. ||, where we have again used reduced units. Here- 
after we will use x to refer to the reduced length, and d 
for the reduced bond distance of the wires (see the inset 
of Fig.§). 

First, we will briefly comment on the data for Au, 
displayed in Fig [j]a. They present two peculiar behav- 
iors: i) The planar zigzag is more stable than the lin- 
ear geometry, even before the latter reaches its optimum 
bond length, and it) the zigzag structure presents two 
energy minima as a function of the wire length. While a 
zigzag structure is naturally expected under compression, 



5.29 A, 4.3 GPa and, 0.94 eV/atom for K. 



the first point mentioned above implies that the zigzag 
distortion, which increases the bond length, is favorable 
even when the linear wire is stretched. This is quite sur- 
prising in a metallic system like gold where, in principle, 
we do not expect any important bond directionality ef- 
fects. Fig. [l]b shows that this result also holds for Cu 
wires. As clearly shown in the insets of Fig. |l|a and b, in 
both cases the linear configuration is, at its equilibrium 
length (x=l), already ~0.1 eV/atom less stable than the 
zigzag structure. In spite of this energy gain, Fig. || shows 
that for both materials, at x=l, the nearest-neighbors 
distance increases in ^5 % when going from the linear 
to the zigzag conformation. The behavior of the energy 
shown in panels (c) and (d) of Fig. [|, corresponding to Ca 
and K, is quite different. In those curves we do not find 
any sign (at least within the precision of the calculation) 
of the zigzag distortion until the linear chain starts to be 
compressed. This is what we could call a 'trivial' zigzag 
distortion: the energy minimum of the linear wire ap- 
proximately becomes an inflection point for the energy 
of the zigzag configuration. This is in fact the kind of 
behavior predicted with, for example, a simple pairwise 
potential. 

All the curves in Fig. [l] present a clear minimum 
near x ~ 0.5. This minimum corresponds to a zigzag 
formed by equilateral triangles, where each atom has 
four nearest-neighbors at the same distance. In Fig. ^ 
the equilateral configuration corresponds to the points 
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FIG. 2. Evolution of the bond length d for monatomic 
wires with a zigzag geometry, as a function of the length per 
atom x. All lengths are expressed in units of the equilibrium 
bond length of the linear wires. Open and solid diamonds 
stand for Au and Cu, respectively, circles for K, and squares 
for Ca. The dashed line indicates the distance 2x between 
second- neighbors. 
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FIG. 3. Band structure of an infinite monatomic copper 
wire in a linear geometry of length 2.20 A/atom (a; =0.98) 

(a) , and a planar zigzag geometry of 1.90 A/atom (a; =0.84) 

(b) , and 1.30 A/atom (x =0.58) (c). 



where each curve intersects the dashed line, i.e. where 
d = 2x. 

The second peculiarity of the Au wires mentioned 
above, i.e. the presence of two minima, is not observed 
for the other species studied. In ref. || we explained this 
behavior with a simple jellium model, enphasizing the 
relationship between the energy vs length behavior, and 
the crossing of the Fermi level by different bands. For 
linear wires of Cu and Au with x < 1.1, a doubly degen- 
erated 7r band (with mainly d character) rises through 
the Fermi level (Fig. |]a). In the case of Au linear wires, 
this was already reported in ref. ]3j , and later confirmed 
in ref. 0. The high density of states at the Fermi level, 
associated to this band, destabilizes the linear chain, and 
is the origin of the 'anomalous' zigzag behavior of Au and 
Cu (as opposed to the 'normal' behavior of Ca and K), 
which begins when the crossing occurs. 

For x > 0.7, Au, Cu, and K zigzag-wires have a single 
band crossing of the Fermi level (Fig. ||b) , at the Brillouin 
zone edge. This situation might, in principle, favor a 
Peirls distortion. However, for the wire lenghts studied 
here, we have not found any measurable dimerization. As 
the wire is contracted below x ~ 0.7 (the exact positions 
are indicated by arrows in Fig. ^) a new conduction band 
crosses the Fermi level. This is shown in Fig. ||c for the 
Cu wire. For the three elements (Au, Cu and, K), this 
new band crossing precedes a change in the behavior of 
the energy, and may be interpreted as the begining of a 
chemical bond between second neighbors. 

The case of the Ca zigzag- wire is different. The wire 
is semiconducting for most of the length interval studied 
here, and the bands cross the Fermi level at a much lower 
x ~ 0.45 (below the equilibrim length) making it metal- 



lic. At x ~ 0.75 its band structure suffers another qual- 
itative change: the gap, that was direct at T for larger 
lengths, begins to move towards the edge of the two-atom 
Brillouin zone, still being approximately direct. 

While the energy curves of K and Ca zigzag- wires were 
relatively similar, the behavior of the bond distance, dis- 
played in Fig. ^, is completely different for metallic and 
semiconducting wires. In the case of the metallic wires 
(Au, Cu and, K) the bond distance increases as the zigzag 
wire is contracted. This can be understood as a conse- 
quence of the metallic bonding: as the second-neighbour 
distance decreases, the first-neighbor distance must in- 
crease to maintain the optimum value of the effective 
atomic coordination. In fact, the calculated behavior is 
qualitatively reproduced by semiempirical metallic po- 
tentials like the glue model ]l9| , or the embedded atom 
method [^o| , or even a jellium model where the volume 
per atom remains fixed. However, for the semiconducting 
wire (Ca), the nearest-neighbor distance behaves just in 
the opposite way, decreasing as the wire is compressed. 
It is worth noting here that the predictions of a simple 
pair potential, like the Lennard- Jones pi} ] , lie exactly in 
between these two regimes, i.e. the nearest-neighbors 
distance would remain almost constant once the wire 
assumes the zigzag configuration. In Fig. g, it is also 
clearly visible the abrupt change in the behavior of the 
bond length of K wires induced by the band crossing at 
x ~ 0.7. 

In conclusion, we have applied first-principles density- 
functional calculations to study the energetics, and the 
elastic and electronic properties of monatomic wires of 
different elements, in a linear and a planar zigzag geome- 
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try. Surprisingly, for Cu and Au wires, a zigzag distortion 
is favorable even when the linear wire is stretched. This 
phenomenon is not observed for K and Ca wires, and 
might be related with the presence of d bands (with a 
high density of states associated) at the Fermi level for 
the linear conformation of Au and Cu wires. For all the 
studied materials, the equilibrium structure is an equi- 
lateral zigzag (bond angle of 60°). Only in the case of 
Au, the zigzag geometry can also be stabilized for an 
intermidiate bond angle of 131°. The behavior of the 
nearest-neighbor bond length is qualitatively different in 
the metallic (Au, Cu and, K) and semiconducting wires 
(Ca). 
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